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Abstract 

We investigate the folding transition of a single diblock copolymer consisting of a semiflexible and a 
flexible block. We obtain a Saturn-shaped core-shell conformation in the folded state, in which the flexible 
block forms a core and the semiflexible block wraps around it. We demonstrate two distinctive features of 
the core-shell structures: (i) The kinetics of the folding transition in the copolymer are significantly more 
efficient than those of a semiflexible homopolymer. (ii) The core-shell structure does not depend on the 
transition pathway. 
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I. INTRODUCTION 



Polymer-based nanostructures have been extensively studied due to their importance in indus- 
trial applications, particularly in nanodevices and nanomachine*^. They are also of importance 
in living cells; DNA and protein molecules have nano-ordered structures, which show a close 
relationship with their biological functions. Macromolecules in biological systems typically un- 
dergo conformational transitions. In polymer physics, this has been discussed as the coil-globule 
transition 4 , in which a flexible swollen homopolymer collapses into a spherical globular conforma- 
tion. The globule state is liquid-like and disordered, whereas a lot of macromolecules, particularly 
biomacromolecules, have ordered, folded structures. For this reason, there has recently been much 
attention to semiflexible homopolymers, which show bending rigidity along the chain, and, as a 
result, have a rod-like properties, although their contour lengths are long enough to exhibit overall 
fluctuations. A semiflexible homopolymer has been extensively studied as a model for a DNA 
molecules^^. In these works, nano-ordered structures, such as toroids, cylinders, and rackets 
were investigated experimentally, theoretically and computationally 8 . Due to these extensive stud- 
ies, it was found that the formation of nano-ordered conformations in DNA is well reproduced in 
semiflexible polymers with homogeneous bending rigidity. In living systems, however, biopoly- 
mers such as proteins are, in general, heteropolymers with complicated sequences of amino acids. 
It is often mentioned that sequences within proteins are relevant for their conformations in the 
folded states. Therefore, most of studies heretofore conducted have dealt with sequences in het- 
eropolymers which lead to diversity in conformations. Nevertheless, our understanding of role of 
heterogeneity in bending rigidity is still primitive. In this article, we propose a minimum model to 
extend the concept of semiflexible polymers toward heteropolymers. A single block copolymer is 
the simplest extension of a single homopolymer in the direction of single heteropolymers such as 
proteins. To this end, we concentrate on a simple model: diblock copolymers in which two blocks 
possess different levels of flexibility. 

Assemblies of rod-coil copolymers have been first discussed in the context of polymer 
blends^; they have recently been investigated in the context of polymer solutions 11 , with re- 
gard to their applications in the field of nanocapsules^. However, less attention has been paid 
to the properties of single rod-coil copolymers. In computer simulations, several novel structures 
were found by Cooke et al. for multi-block rod-coil copolymers, which they refer to as semi- 
flexible copolymers. The static properties of these species were studied 13 . In the present study, 
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we focus instead on the kinetics of the process of folding into ordered phases in single rod-coil 
diblock copolymers. 

The kinetics and pathways of conformational transitions are relevant since it is well known 
that in experiments the results strongly depend on the method of preparation. For example, the 
authors demonstrated that the final structure in the folding transition of DNA depends on whether 
condensing agents or monovalent salts are added first= This result suggests that we must discuss 
not only static stability but also kinetics pathways for the folding transition. Along these lines, we 
have investigated the dependence of folded structures of semiflexible polymers depend on kinetic 
pathways^. Once a chain falls into a metastable state, it cannot escape to the most stable state 
under thermal fluctuations. For semiflexible polymers, even when the toroidal conformation is the 
most favorable, a small hairpin-like conformation in the early stage of the folding transition results 
in the formation of a cylindrical structure as a stable state. Due to a high free-energy barrier, the 
cylinder does not make the transition to a toroid over a practical time scale. Our interest is in 
the kinetics of formation of specific nano-ordered structures, and particularly on the pathways to 
folded states. We demonstrate that our single flexible-semiflexible copolymer avoids local free 
energy minima and forms a core-shell structure. 



II. SIMULATIONS 



We carried out Langevin dynamics simulations for a single rod-coil block copolymer. A rod- 
coil copolymer consists of a flexible and a semiflexible block. Although rigid rods with infinitely 
large bending elasticity are often used for the rod block, in order to allow more general discussions, 
we used a semiflexible polymer with bending rigidity. We adapted a bead-spring model with the 
following potentials: 



Vbeads = 2 / A \\ T i+i ~ V i\ ~ a ) ' 

i 

Vbend = 2^K a 2](l -COSdf), 
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where r, is the coordinate of the zth monomer and 6i is the angle between adjacent bond vectors. 
The subscripts a and fi denote "s" (semiflexible) or "f" (flexible) monomers. The monomer size 
a was chosen as the unit length, and k B T as the unit energy. The excluded volume and short- 



ranged attractive interactions between monomers are included with the Lennard- Jones potential 
with the coupling constant e a p, which determines the strength of the attractive interaction between 
monomers in the a and /S states. We set the spring constant to be k = 400. The bending elasticity 
for a flexible block was chosen to be *y = 0. We note that the persistence length of a semiflexible 
block can be written as l p T = K s a. The length of each block was chosen as N s = Nf = 128. 
The equation of motion is written as, 

d 2 r ; - dr, dV 

where m, y are the mass and friction constant of monomeric units, respectively. The unit time 
scale is r = ya 2 /k B T. We set the time step as O.Olr, and use m = 1.0 and y = 1.0. With these 
parameters, the relaxation time of the momentum of a monomer is sufficiently fast compared to 
the time scale of interest. Gaussian white noise £ satisfies a fluctuation-dissipation relation, 

< £.(f) • £..(f') >= 6yk B T6 u 6(t - f'). (5) 



III. CORE-SHELL STRUCTURE 



First, we consider the system in non-selective solvents, i.e. e ss = e s f = e/y = e. Figure 1 
shows typical conformations of folded single rod-coil block copolymers. At sufficiently large e, 
a rod-coil copolymer in the elongated state undergoes a transition to the folded state. For small 
value of k s , the collapse is disordered, while a Saturn-shape conformation is obtained when k s is 
large. In this state, the flexible block at the core is surrounded by the semiflexible block (Fig. Q] 
(II)). Intuitively, our core-shell structure is a composite of a globule and toroid. This is a typical 
characteristics of block copolymers: the folded state of this block copolymer incorporates the 
characters of both flexible and semiflexible polymers. This conformation is observed over a wide 
range of parameters, while at k s » and e <sc 1, the semiflexible block is unable to fold, resulting 
in a tap-pole conformation with a long lifetime (Fig. ED (III)). 

The core-shell structure is robust against the change of the length of blocks. The flexible block 
forms spherical shape, and thus, the size Rf is proportional to 

R f ~N 1 / 3 . (6) 

Since the exponent is much smaller than 1, the size is insensitive to the length. Furthermore, the 
size of the semiflexible block, R s is 

^~^ 1/5 , (7) 
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which is also insensitive to the length of the block. In fact, both Nf = 64 and Nf = 256 blocks 
reproduce the core- shell structure in simulations (data not shown). 

The final conformation of the rod-coil copolymer is the uniquely determined core-shell struc- 
ture, in which the semiflexible block can only assume a toroidal conformation. In contrast, a 
semiflexible homopolymer possesses metastable conformations, such as a cylinder (fig. HP). The 
toroidal and the cylindrical conformations are both stable, and, over the practical time scale, the 
distribution does not change since a transition between two states does not occur. We will come 
back to this point later. 

We show in Fig. 2 the kinetics of the formation for a Saturn-shaped structure, where we sud- 
denly increase e at t = 0. The flexible block collapses at an early stage, and the semiflexible 
block then gradually wraps around it. This is contrast to the situation for flexible and semiflex- 
ible homopolymers: flexible polymers exhibit spinodal decomposition in the folding process^, 
whereas semiflexible polymers undergo a nucleation and coarsening process-. Due to a stochas- 
tic feature of nucleation, there is a long lag time before the folding transition in a semiflexible 
polymer takes place. These observations indicates that a flexible polymer collapses much more 
quickly than a semiflexible polymer. Figure 2 shows the time evolution in the collapse ratio. As 
we can see, rod-coil copolymers undergo folding more quickly than semiflexible homopolymers. 
While at t ~ 350, a rod-coil copolymer undergoes complete folding on average, only 30 % of the 
monomers in a semiflexible homopolymer exhibit a collapsed states. The folding of semiflexible 
homopolymers is subject to a long lag time originating from the nucleation process, while rod-coil 
copolymers start to fold without such a process due to the existence of the flexible-block core. As 
a result, the folding transition proceeds quickly. 

These results are confirmed by the consideration of the time scale of folding in flexible and 
semiflexible polymers 17 . The initial nucleation time depends exponentially on the persistence 
length of the polymer. This leads to a clear separation of the time scale between the flexible and 
semiflexible blocks. Thus, the folding kinetics of a rod-coil copolymer can be expected to be a two- 
step process: folding of the flexible block, and coarsening due to adsorption of the semiflexible 
block onto the globule of the flexible block. 
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IV. PATHWAYS 



Next, we move our focus to pathways to the folded state. Although in the above case all 
the parameters were quenched simultaneously, this is not always the case. For this reason, we 
consider two-step procedures for changing the parameters, {e ss , e S f, ejf} from {0.30,0.30,0.30} 
to {1.0, 1.0, 1.0}. There are six possible intermediate states, as shown in Fig. [31 which have 
various intermediate conformations depending on the procedures used. Several of the interme- 
diate states show cylindrical shapes in their semiflexible parts. Therefore these structures could 
be considered to be metastable states because the core-shell shape mentioned earlier contains a 
toroidal conformation. However, this is not the case. Interestingly, the final structure does not 
depend on the procedures used. This indicates that our core-shell structure is also independent 
of the order in which the procedure takes place. This point is clearly demonstorated in the case 
of {e ss , e s f, 6ff} = {1.0, 1.0,0.30} (Fig. [3F), where AR = R s - R f changes negative to positive. 
Under these conditions, the intermediate state exhibits a core-shell structure in which the core is 
the semiflexible block and the shell is the flexible block. As we showed, in the final core-shell 
structure (Fig. 1(H)) the core is the flexible block and the shell is semiflexible block; therefore, the 
intermediate state shown in Fig. 3(F) has an inside-out composition. Nevertheless, the core-shell 
shown in Fig. 1 (II) is always obtained as the final structure. 

Why is the core-shell conformation independent of the pathway? First, we should stress that a 
semiflexible polymer has several metastable states, including a toroid and a cylinder. As shown in 
Fig. [Ip, the conformation of a semiflexible homopolymer depends on pathway of the transition. 
The formation of a cylindrical shape at an early stage of the kinetics leads to a cylinder at the final 
state. The free energy barrier between the states is sufficiently high that the transition between 
them does not occur over a practical time scale. However the addition of a flexible block removes 
the metastable states. Hereafter, we discuss the reasons for pathway-independence. To see the 
stability of the cylindrical conformation, it is reasonable to assume the size of the cylinder is of the 
order of the persistence length, because cylinders obtained experimentally and in simulations are 
around this size. In this case, the conformational fluctuation is too small to undergo a transition 
to a toroidal shape when two ends of a cylinder meet. The only pathway from a cylinder to a 
toroid is to make a hole at the center of a cylinder. We consider a complex of a globule composed 
of a flexible block and a cylinder composed a semiflexible block. When the globule is slightly 
embedded in the cylinder (Fig. |4}, the free energy of the system due to the deformation, compared 
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with the reference state, is 

AF = F mt + F sur + Fbend- (8) 

The reference state is chosen as the conformation in which the globule and the cylinder are sep- 
arated, without any interactions. F int , F sm and Fbend represent the attractive interaction between 
the flexible and semiflexible blocks, and the increase in surface energy and the bending energy 
due to deformation. The energy represented by the first term contributes to deformation, while 
those of the second and third terms prevent it. Since the attractive interaction is short-range, F int 
is proportional to the contact area. This energy per unit area can be represented by -e S f/a 2 , and 
thus, the free energy is written as 

Ane sf R 2 f 

Fm r^ Ae /) 2 ' ( 9 ) 

a 1 

where A8f is the angle characterizing the deformation (Fig. IJ). The bending free energy at the 
lowest order in A6 f is proportional to square of the curvature, which is approximately A8 f /L, and 
thus 

Fbend " oi(^r) (A0/)4, (10) 

where L is the length of the longer axis of the cylinder. With volume conservation and the condition 
Rf sin 8f = R s sin 9 S , we find that the surface free energy is proportional to (A9/) 4 : 

F sur - ^(A#/) 4 - (11) 

The fact that F int is dominant for small A6f indicates instability in the cylindrical conformation 
and an increase in contact area. We thus obtain the core-shell conformation as the equilibrium 
state of the global free energy minimum. 

V. SUMMARY AND REMARKS 

We report novel nanostructures made from single rod-coil block copolymers. The core-shell 
structure is found to be the equilibrium state. This structure is obtained due to incommensuration 
of the flexible and semiflexible blocks: the segments in semiflexible blocks tend to align in folded 
states, while those in flexible blocks tend to be disordered. This state is the counterpart of the 
coexistence states found in rod-coil copolymer blends^. 

The kinetics of the block copolymers are markedly different from those of a semiflexible ho- 
mopolymers: in a block copolymer, the flexible blocks tend to collapse quickly, which means that 
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the semiflexible blocks can collapse without a nucleation process by making use of the globular 
collapsed parts of the flexible blocks. On the other hand, semiflexible homopolymers require a 
long time to enter the folded state because the time required for nucleation appears as a lag-time 
before the collapse. We also discussed the robustness of the core-shell structure: formation of the 
core-shell structure does not depend on intermediate structures. This was discussed in terms of the 
instability of cylindrical shapes. 

Although our simple model does not specify real macromolecules to which it can be applied, 
inhomogeneity of elasticity is realized, for example, with binding ligands^, and in mixtures of 
a helix and a coil 2 ^ where the helix is relatively stiff. In protein folding, it was proposed that 
secondary structures are transiently formed at an early stage according to their local preference, 
independently of native structures 2 ^' 3 .. This leads to inhomogeneity of semiflexibility along a 
chain. Thus, we expect that general features of kinetics at later stage of proteins in this class 
are to be described with a mesoscopic model of inhomogeneous bending rigidity. In fact, our 
polymer has unique folded structure regardless of difference in pathways. Mesoscopic models for 
protein folding are still developing. It was proposed that correlation in sequences of hydrophobic- 
ity is relevant to the core-shell structure and, in fact, they reproduced the structure with designed 
sequences 24 . Our approach focuses instead on semiflexibility along a chain, and we expect that 
inhomogeneity in semiflexibility plays a key role for transition kinetics in biomacromolecules. 
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List of Figures 

1 Phase diagram of single folded rod-coil copolymers at various e and k s (A), and 
illustrations of the conformations the polymer in each phase: (I) disordered col- 
lapse, (II) Saturn-shape phase separation, (III) tadpole, and (IV) swollen. The 
probability distribution of the gyration radius of a rod-coil copolymer (B) and 
a flexible (C) and semiflexible (D) homopolymer in the collapsed state are also 
shown. In the figures, the monomer in the flexible block is shown in light gray, 
while the monomer in the semiflexible block is shown in dark gray. A semiflexible 
homopolymer has two stable conformations, toroidal and cylindrical 

2 Time evolution of the number of collapsed monomers in the semiflexible block of 
a rod-coil copolymer (O) and a semiflexible homopolymer (A). Each polymer has 
N = 256 monomers, and a bending elasticity k s = 15 for the semiflexible parts. 
We change the attractive interaction at t = to e = 0.8. iVfoid shows the nuber 
of monomers in the folded state. The bare data are shown in (A), while the mean 
evolution over 30 runs is shown in (B). Snapshots of the polymer conformation 
during the transition for a rod-coil copolymer (a-f) and a semiflexible homopoly- 
mer (a'-f ') are shown in (C) 

3 Pathways of the folding kinetics of a rod-coil copolymer. The effective sizes 
of the flexible (a = f) and a semiflexible (a = s) blocks, R a are defined as 
R a = J rp a (r)dr/ J p a (r)dr, where p a is number density of each block in the radial 
direction. Here, r represents the distance from the center of mass. The value of 
AiR = R s - Rf, which represents the extent to which the flexible block is located 
outside the semiflexible block, is plotted for the various intermediate states. The 
systems are quenched into six intermediate conditions at t = 0, and then into the 
final condition at t = 2 x 10 6 as indicated by the dashed line. Representations 
of the initial, intermediate and final configurations are shown in the figure. The 
trajectories are averaged over 100 runs 

4 Deformation of a cylinder resulting from interaction with a globule composed of a 
flexible block. The left-hand figure shows a cross-section of the complex consist- 
ing of a spherical globule and a cylinder 



10 



Figures 



A. 



I. • 

0.6 
■ ■. I 



{ID • 



i . inn 



TV 

.A 



5.0 10.0 15.0 200 



• I! 



.. n 

II 111 



inn 



3 ID U Zl D J ID Zi 31 D J ID 25 Z] ZS 



FIG. 1: Phase diagram of single folded rod-coil copolymers at various e and k s (A), and illustrations of the 
conformations the polymer in each phase: (I) disordered collapse, (II) Saturn-shape phase separation, (III) 
tadpole, and (IV) swollen. The probability distribution of the gyration radius of a rod-coil copolymer (B) 
and a flexible (C) and semiflexible (D) homopolymer in the collapsed state are also shown. In the figures, 
the monomer in the flexible block is shown in light gray, while the monomer in the semiflexible block is 
shown in dark gray. A semiflexible homopolymer has two stable conformations, toroidal and cylindrical. 
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FIG. 2: Time evolution of the number of collapsed monomers in the semiflexible block of a rod-coil copoly- 
mer (O) and a semiflexible homopolymer (A). Each polymer has N = 256 monomers, and a bending elas- 
ticity k s = 15 for the semiflexible parts. We change the attractive interaction at t = to e = 0.8. A^foid 
shows the nuber of monomers in the folded state. The bare data are shown in (A), while the mean evolution 
over 30 runs is shown in (B). Snapshots of the polymer conformation during the transition for a rod-coil 
copolymer (a-f) and a semiflexible homopolymer (a'-f ') are shown in (C). 
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FIG. 3: Pathways of the folding kinetics of a rod-coil copolymer. The effective sizes of the flexible (a = f) 
and a semiflexible {a - s) blocks, R a are defined as R a - J rp a (r)dr/ J p a (r)dr, where p a is number 
density of each block in the radial direction. Here, r represents the distance from the center of mass. 
The value of = R s - Rf, which represents the extent to which the flexible block is located outside 
the semiflexible block, is plotted for the various intermediate states. The systems are quenched into six 
intermediate conditions at t = 0, and then into the final condition at t = 2 x 10 6 as indicated by the dashed 
line. Representations of the initial, intermediate and final configurations are shown in the figure. The 
trajectories are averaged over 100 runs. 
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FIG. 4: Deformation of a cylinder resulting from interaction with a globule composed of a flexible block. 
The left-hand figure shows a cross-section of the complex consisting of a spherical globule and a cylinder. 
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